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^ Abstract 
< 

Our interest is in a class of directed solid-on-solid models, which may be regarded as 
p<| continuum versions of boxed plane partitions. In the case that the heights are chosen 

from a uniform distribution, the joint PDF of the heights is the same as that for the 
positions in a unitized bead process recently introduced by the authors and Nordenstam. 
We use knowledge of the correlation functions for the latter to show that upon a certain 
scaling the fluctuations of the heights along the back row of the solid-on-solid model are 
given by the Airy process from random matrix theory, as is the case for boxed plane 
partitions. Moreover, we show that this limiting distribution remains true if instead of 
the uniform distribution, the heights are sampled from a general absolutely continuous 
distribution. 

> 

1 Introduction 

in 

For more than a decade now, it has been demonstrated by numerous examples that fluctuation 
formulas relating to the largest eigenvalues of certain random matrix ensembles are also charac- 
teristic of fluctuation formulas for certain random growth models (for a recent review see [7j). 
Because the precise functional form of the random matrix quantities are known (see e.g. [9|), 
there are thus precise theoretical predictions for the functional form of the corresponding quan- 
tities in the growth models. Very recently, the validity of these theoretical predictions [TB] has 
been demonstrated by a high precision experiment involving curved interface fluctuations of a 
growing droplet [T9] . 

Crucial to the theoretical predictions is the hypothesis of universality of the fluctuation 
formulas: they depend on some global symmetry properties of the interface (in particular, is it 
curved or flat?), but not the details of the microscopic interactions. This hypothesis is generally 
unproven, as the results for many of the random growth models rely on the integrability of the 
particular microscopic choice of parameters. One recent advance has been to show that the 
random matrix fluctuation formulas predicted on the basis of exact calculation have been shown 
to hold true for the continuum KPZ equation with narrow wedge initial conditions [U [T7j. But 
mostly universality results are out of reach. 
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An outstanding problem of this type relates to last passage percolation. Let {<%} be 
independent random variables drawn from a probability distribution with unit mean and finite 
variance. The random variable 

In = max a^- (1-1) 

(l,l)u/r(jV,AT) 

finishing at (N,N) and forming an up/right directed path when plotted on the integer grid, 
is referred to as the last passage time. In the case that each is drawn from an exponential 
distribution (or more generally a geometric distribution) the scaled fluctuation of In about its 
mean AN is known to be identical to the fluctuation of the scaled largest eigenvalue of a random 
complex Hermitian matrix [T3] . But for the with distributions different from these solvable 
cases, the corresponding theorem is not known. 

In this paper we will demonstrate universality for fluctuations of a certain two-dimensional 
directed solid-on-solid model, defined in terms of a grid of height variables [:%]y=i,...,Ar- Each 
height variable is chosen independently from the same non-negative, absolutely continuous 
distribution P, subject to a constraint on its value relative to its neighbours. The solid-on- 
solid model is defined in Section 2. In Section 3, in the particular case T> = Umi] (the uniform 
distribution on [0, 1]), use is made of a correspondence between this instance of the solid-on-solid 
model and the recently introduced finitized bead process [8], to obtain a Fredholm determinant 
formula for the distribution of the largest heights along the back row of the former. In Section 
4 the explicit form of the distribution function between scaled heights on the back row 0(iV 2 / 3 ) 
away from each other is computed. The functional form found is precisely that for the Airy 
process [TJ)J [TTJ [TB], which in turn is identical to the correlation functions for the GUE minor 
process at the soft edge \±0\. We show in Section 5 that, upon an appropriate choice of scaled 
variables, this same correlation function persists for the Xij chosen from a general absolutely 
continuous distribution T>, thereby establishing the claimed universality. 



2 Definition of the model 

Consider the N x N integer grid : 1 < i,j < N}. At each site associate a 

height variable Xi.j = T>, where D is an absolutely continuous distribution with density h(x). 
Furthermore, require that for fixed % 



Xi,l < %i,2 < ■ ■ < %i,N 

(heights increase along rows) and that for fixed j 

X\j <C %2,j ^ ' ' ' ^ XNJ 



(heights increase up columns). Thus a directed up/right lattice path (recall text below (1.1)) 
must encounter successively larger heights. 

Rotate the square grid by 45° anti-clockwise and mark lines parallel to the y-axis through 
the lattice points. For aniVxiV grid there are 2N — 1 lines, and on the l-th line there are I 
grid points for I — 1, . . . , N, and 2N — I grid point for I = iV + 1, . . . , 2N — 1. Let yf denote 
the height of the i-th lattice point counting from the top on line j. We have that 1 < i < N(l), 
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where N(l) := I, (i = 1, . . . , N), and N(l) := 2N - / (Z = N + 1, . . . , 2N) in keeping with the 
total number of particles on line j, the relationship 



(N+l-j+k) . , 

y(N+i-j+k)+i-k-> j — 

(N+l-j+k) • , 

»(;v+i-i+fc)+i-i' J i 



and most importantly the interlacings 



and 



< yr i> < Vj', l<j<2N-(l + l)<{N- 1). 



(0 / C-i) ^ (0 

y)U<y) } <yy, 



i < i < (i - 1) < iv 



(2.1) 



(2.2) 



Let us denote the interlacings (2.1) and (2.2) by Y({y}). The joint PDF for a configuration 



specified by the height coordinates ^f=i{Vj }j=i,...,N(i) will then be give in terms of the PDF 



for the heights h(x) by 

^ 2N N(l) 
1=1 j=l 

where xa — 1 f° r A true and \a = otherwise. 



(2.3) 



As our first exact result in relation to the joint PDF (2.3), let us compute the normalization 



C. 



Proposition 2.1 The normalization C in (2.3) is independent of the PDF h, and is given by 



N 



r(i+£ 



3=0 



r(N+j + iy 



(2.4) 



Proof. The PDF h(x) maps lR>o to K>o- Its antiderivative 



H(x) = / h(t) dt 
Jo 

is therefore a non-decreasing function from IR>o to [0, lj.Thus if we change variables 

uf = H{yf) (2.5) 
we see that the interlacings (2.1), (2.2) remain true in the variables {uf^}- Furthermore 



du 



h = Kvf) 



a) 



and so 



dy 



2N N(l) 2N N(l) 

n n Kvf)xr{{ V }) dy= n n u m {uf)xv{{u}) du. 

1=1 3=1 1=1 3=1 



(2.6) 
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The measure on the RHS of (2.6) is an example (the case p = q = N) of the unitization 
of the bead process introduced in [5]. But note that the interpretations of the variables are 



different: in [8 J Uj is the position of the particle j on line I while here u y - } is a (transformed) 
height at a lattice point as specified by (2.5). 



(i) 



Also given in [8] is a discrete version of (2.6), and furthermore summation methods are 
presented to compute marginal PDFs for a single line. In [9, Prop. 10.2.3] this working is 
adopted to the present continuous setting. In particular, it is shown (HI Eq. (10.66)] that 



integrating over lines I = 1,2, ... , with configuration {yj +1 ^}j=i t ...,N+i gives 



N 1 

n 

s=1 l<j<k<N+l 



(y 



(N+l) 



y f r ih , 



By symmetry, this same expression results by integrating over lines / = 2N, 2N 



(in that order) with configuration {yj N+1 ^} 



j=l,...,N+ 



i on line N + l. Hence 



C 



N 

n 

8=1 



1\ 2 



n <* 



l<j<k<N+l 



(2.7) 
N + l 

(2.8) 



where Y denotes the region 



l>2/i 



(N+l) 



> 2/2 



(N+l) 



>-">Vn 



(N+l) 



> 0. 



The integrand in (2.8) is symmetric, so the region Y can be replaced by [0, 1] N+1 provided 
we divide by (N + 1)!. Furthermore, 



[0,1] 



JV+l 



n <* 



N 



(N+l) 



Vk 



(N+l)^2 



dj/i»+» = n 



l<j<k<N+l 



3=0 



r 2 (i + J )r(2 + j; 

T(N + j + l) 



(this is a special case of the Selberg integral (see e.g. (HJ Eq. (4.3)]), or alternatively can be de- 
rived using knowledge of the normalization of single variable Jacobi polynomials [9, Eq. (5.74)]) 
and (Q results. □ 



3 Uniformly distributed height variables 

Suppose the height variables are sampled from the uniform distribution on [0,1]. It was re- 



marked during the proof of Proposition 2.1 that another interpretation of the joint PDF (2.3) 
for the heights in the directed solid-on-solid model is then as a joint PDF for interlaced par- 
ticles on 2N — 1 parallel line segments [0, 1], with N(l) particles on line /. It was furthermore 
remarked that this latter model is the so called finitized bead process, studied in [8], or more 
precisely the special case 

p = N, q = N (3.1) 

of the finitized bead process (general p and q would correspond to the height model of Section 
2 being defined in a rectangle rather than square). 
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Known results for the bead model can then be used to deduce properties of the directed 
solid-on- solid model. Of particular interest for present purposes are the facts that configurations 
of the bead process can be sampled by generating eigenvalues of a random matrix process, and 
that the correlation functions have an explicit determinantal form. We will revise these results 
separately, and discuss their consequence in relation to the directed solid-on-solid model. 



3.1 Sampling using random matrices 

What does a typical configuration of the directed solid-on-solid model look like? Section 4 of 



[H] details how (a generalization) of the PDF (2.3) in the case that the heights are sampled 



from the uniform distribution C/[o,i] can be obtained as the joint eigenvalue PDF of a sequence 
of random matrices. For a given matrix A, the construction is based on random corank 1 
projections, M = UAH. Here II = I — fr , where x is a normalized complex Gaussian vector of 
the same number of rows as A. The effect of the random corank 1 projection is to reduce by one 
the multiplicity of any degenerate eigenvalues in A, and also to create a new zero eigenvalue. 

The first step is to form Mi = nAiIl, where A\ = diag((0) 7V , (1)^) (the notation (k) is 
used to denote k repeated N times). The eigenvalue of Mi then corresponds to the height h$.. 

Next, for r = 2, . . . , N, inductively generate {/i^ }i=i,...,r as the eigenvalues different from and 
1 of the matrix M r := IL4j.II, where 



A r = dmg((0) N - r+ \ht 1 \...,hti\(l) 



N-r+l 



And after this, for r = 1, . . . , N — 1, inductively generate {/i^ - }i=i,...,jv-r (r = 1, . . . , iV — 1) 
as the eigenvalues different from of M N+r = UA N+r Il where 



A N+r = di&g(h[ N+r 1 \...,h'^ +r 1} 
Crucial to the practical implementation of this construction is the fact that if 

A = diag((ai) Sl ,(a 2 ) 



the n — 1 eigenvalues different from {aj}j = i r .. in and occur at the zeros of the random rational 
function 

n 

Qi 



E 



x — a, 
i * 



where (qi, . . . ,q n ) has the Dirichlet distribution D[si, . . . , s n ] [T2] . 

We have made use of the above theory to generate a typical configuration of the directed 



solid-on-solid model in the case p = q = 10, which is displayed graphically in Figure 3.1 
We remark that if instead of sampling the individual heights from the continuous uniform 
distribution, they were sampled instead from the discrete uniform distribution on [0, c] say 
(c G Z + ), then our directed solid-on-solid model would be equivalent to boxed plane partitions, 
of box size N x N x c. There is no longer a random matrix approach to the sampling of a 
typical configuration, but we point out that one recently introduced method jl] (see also [2] 
and |5J) has the interpretation as defining a dynamics on the underlying interlaced particle 
configurations. 



5 



Figure 1: A typical configuration of the directed solid-on-solid model on a 10 x 10 grid, with the 
heights scaled by a factor of 10. This is constructed using eigenvalues of a particular sequence 
of random matrices as described in the text. 



3.2 Correlations and distribution functions 

Let (/, x) refer to the position x on line I of the bead process. The correlation function p^(X) 

for the configuration X = Uf =1 {Zj, . . . , xi\ *'}, where each x^ is on line U and k = Y^i=i r^, 
defined as the average 

pwpo=(n e n^-vS)) (3.2) 

i=l mi/-^tn ri <J=1 

and thus requiring that k of the particles be fixed at positions X, has been computed in [8j 
Prop. 5.1]. The correlations are given in terms of the rescaled Jacobi polynomials 

P^ b \x):=Pi a > b \l-2x), (3.3) 

which satisfy the orthogonality 

x a (l - x) h P^ b \x)P^ b \x) dx = A/f'% fc (3.4) 



where 



u{a , b ) = 1 (n + a)\(n + b)\ 

2n + a + b+l n\(n + a + b)\ ' 1 ' 
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We also require the quantities 



a ^ = {l, N + l<l<2N-l (3 - 6) 

/ (-1)*"', 1 < / < N 

bl{X >- \ x l - N (l-x) l - N , N + l<l<2N-l {6J} 

'* (N-j)\ 



Cf = \ J) ; 1 < l < 2N - 1 (3.8) 



Q^-l^^ (3 9) 

^(x)4^— ^ N (3 10) 

(X) " \ Af^~ N \ N + 1 < I < 2N — 1 (3 - 10) 

a(Z, f ) = min{2iV + 1 — 1, I', iV}. (3.11) 
In terms of this notation Proposition 5.1 of [8 J gives 

p( k ){X) = det[K(li,Xi;lj,Xj)) ijj=lt _.. )k (3.12) 

where 



K(s, u; t, v) 



(3.13) 



To make use of knowledge of this bead process correlation function in the context of the solid- 
on-solid model, we must use it to specify the distribution of the position of specific particles. In 
fact it is sufficient to specify E Q ({li; (v,i, l)}t=i n ) — the probability there are no particles in 
the interval (ttj, 1) of line li (i — 1, ... ,n). Thus the PDF p max ({£j, Wi}i=i ., n ) f° r the position 
of the particles with maximum displacements on lines li, i = 1, . . . , n, then follows by partial 
differentiation according to 

P m " x ({/i,«i}<=i,...,») = dui ... du Eo({h;(ui,l)h= K ..,n). (3.14) 

The crucial point from the viewpoint of the corresponding height model is that this PDF is 
identical to the PDF for the maximum heights on the same lines i = 1, . . . ,n. 

We are thus faced with the task of expressing E ({li] (ui, 1) } i=1 in terms of correlation 
functions, which is in fact a standard exercise P §8.1]. First note that according to the 
definition 

n N(k) 

£ »(ft;(^i)}H ^(IHKi-VeM: 

1=1 1=1 
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Expanding the double product and recalling the definition (3.2) of the correlation function then 
shows 



N(h) N(l n ) 



E ({l h (u h l)i=L,...,„) = ^2 "' 



-1) 



miH hTin 



X 



mi=0 m n =0 
1 /-I 

Ci) . . . 



ui 



m\\ ■ ■ ■ m n \ 

-i 



dyf l] ...l dy^-.-l dy™... ^fc ) p CEairai )(ua=iU^ 1 {(i il yf ) )}),(3.15) 



where the term mi = ■ ■ ■ = m n = is taken to equal unity. 

This is a general formula valid for any two-dimensional particle system confined to 2N — 1 
parallel lines, with N(l) particles on line I (I — 1, ... , 2N — 1). It is furthermore the case that 



when the fc-point (k = Y17=i m «) correlation function p^(Y) has a determinantal form (3.12) 



the multiple sum (3.15) can be summed [9j §9.1]. This can be done by defining the nxn matrix 
Fredholm integral operator K{{{li\ (ui, l))}i=i,..., n ) with kernel 



K(x, y; {h; (m, l)}<=i,„., n ) = Xxe(u u i)K(x, k; y, lj)x y e(u 3 ,i) 
We then have that 

E ({k; (Ui, l)}i=i,..., n ) = det ^1 - K({k; {u u l)}i=i,..., n )J 



i,j=l,-,n 



(3.16) 



(3.17) 



where the meaning of the determinant can be taken as the product over the eigenvalues of the 
operator. 



4 Scaled limits for uniformly distributed heights 



4.1 The limiting shape 

In the large N limit the global density on line of the finitized bead process has been computed 



in [Hj. With the parameters as in 3.1, and the label j of each line scaled j i— > j/N =: S and 
thus < S < 2, the support of the density was found to be the interval 

rl 

2 



ly/S(2-S),± + ±y/S(2-S) 



--:[c s ,d s ]. (4.1) 
And after dividing by N, the explicit functional form of the density was shown to equal 



p(i)(y, s) 



i V(ds-y)(y-c s ) 
7T 2/(1 — 2/) 



(4.2) 



To interpret these results in terms of the directed solid-on-solid model with heights sampled 
from U[o t i], we first agree to scale the N x N integer grid by 1/N so that it is an iV x iV grid 
within the unit square [0,1] 2 . It then follows immediately from (4.1) that the limiting height 



S 



profiles, h(x, y) say, along y — 1, x — 1, y — and x = respectively are 



h(x,l) = -(l + y / x(2-x)) 

Kl,y) 
h(x,0) 
h(0,y) 



2 
1 

1 

2 
1 



1 



(4.3) 



We remark that (4.3) exhibits the general symmetry h(x,y) = h(y,x), which in turn follows 
from the symmetry with respect to the x and y directions of the rule for the interlacing of 
heights in the definition of the model. 

To specify the heights at other positions of the unit square we need to work in the coordinates 
which relate directly to those used for the scaled bead process. The variable S used in (4.1), 
for < S < 1, then identifies the line segment in [0, l] 2 starting at (0, 1 — 5*) and finishing at 
(S, 1). Let t, < t < 1 denote the scaled position along this line segment so that 



{x,y) = {tS,l + S{t-l)). 



(4.4) 



We seek the height h(S,t) at position t along line S. A little thought (cf. [6l eq. (3.5)]) 
shows that h(S,t) is characterized by the equation 



ts 



h(S,t) 



u, s) du. 



(4.5) 



Thus in the bead model picture, the RHS of ( 4.5[ ) gives the expected number of particles 
from the start of the line segment, to position h(S,t) along the segment. The LHS says this 
number of particles is equal to tS. Hence particle number tS along this segment is expected 
to be at position h(S, t) on the line. In the solid-on-solid picture, the expected position of a 
particular numbered particle is the expected height at the point in the square corresponding 
to the numbering of the particle. 

We note from (4.2) and (4.1) that pm(y, S) is symmetrical about y — 1/2. It follows that 



1/2 



u, s) du 



ds 



and consequently 



h(S, 



or equivalently in terms of xy-coordinates, 

h(x, 1 — x) 



i 

2' 

1 

" 2' 



'u, s) du 



S 
2' 



< S < 1 



< S < 1 



< x < 1. 



(4.6) 



This supplements the exact profiles (4.3). 
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Figure 2: A plot of the limiting surface of the directed solid-on- solid model as calculated using 
and (ggh. 



Along other line segments of the unit square, we have to make do with (4.5), although the 
integral can be evaluated explicitly. Thus recalling (4.1) and (4.2), use of computer algebra 
allows us to conclude 



ts 



arcsm v 



i — / J i - d 2 

/2 + VI - d 2 (arctan . 

^ V 1 — 



tt/2) 



(4.7) 



where 



w := (2/i(S, t) - l)/d, d = y/S{2-S). (4.8) 

For practical determination of the height profile, we thus choose (x, y) within the unit square, 
then determine S and t according to 



S = 1 + x — y, 



t = x/(l + x — y) 



(4.9) 



as implied by (4.4). Substituting these values in (4.7) we obtain an equation for h, which is 



solved using a root finding routine. A graph of the resulting shape is plotted in Figure 4.1 



4.2 Scaled correlation of heights along y = 1 

The expected value of the heights along y = 1 in the limiting solid-on-solid model is the first 
entry in (4.3). We seek a scaled limit of the joint distribution function (3.14) for k heights 
along this line, which according to (3.17) requires computing the large form of (3.16). In the 
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limiting procedure x and y must first undergo a linear change of scale. Thus after introducing 
the scaled line number Sj, for Xj corresponding to this line we must write 



Xj = XoiSj) +XjaN- 2/3 



(4.10) 



where x (S) := 1(1 + ^/5(2 - 5)), a is independent of N and can be chosen at our convenience, 
and the factor N~ 2 / 3 is chosen so that in the variable X the spacing between the large heights 
on line S is 0(1). We must also choose the spacing between scaled lines to scale with N. Thus 
we write 

Sj = S* + tjrN- 1 ^ (4.11) 
where S* is the limiting line on the scale of a division by N, and ts* is independent of N to be 



chosen for convenience. The motivation for the choice (4.11 ) is that the scaled correlations are 
the 0(1) in the variables {i,}. 



Since from £3.2 the correlations are given in terms of Jacobi polynomials (3.3), we require a 



uniform asymptotic expansion appropriate for the scaling (4.10). The relevant such expansion 



has recently been given by Johnstone [Hj: it applies when the Jacobi polynomial parameters 
a, b increase with n, and when x is centred about the largest (or by replacing the roles of a, b, 
the smallest) zero, with a further scaling chosen so that the spacing between zeros is of order 
unity. For our purposes, due to the relation (3.3), we formulate the results of [H] about the 
smallest zero of Pn a ' b \x), and thus the largest zero of Pi a ' b \x). 



Proposition 4.1 Define variables n n , ip and 7 by 



a — b 



cos 7 



K n = 2n + a + b + 1, cos ijj 
and using these variables define 

f , ^ 3 2sin 4 (^ + 7 ) 

M n = cosV + 7, °l= 2 ■ ■ ■ 

k„ sin ip sin 7 

Then, for x = M n — a n X , one has the uniform asymptotic expansion 



a + b 



(4.12) 



(4.13) 



P (a, b){x) 



- x) a (l + x) b (l - X 2 ) 



Ai(X) + 0(X 



-2/31 



e-*/ 2 , X>0 
1, X < 



(4.14) 



Equivalenty, in terms of the rescaled Jacobi polynomials (3.3), for x — (1 — M n )/2 + a n X/2, 
Pi a ' b) (x) -- 



2x a + 1 (l-x) b + 1 



Ai(X) + 0(N- 2/3 ) 



e- x / 2 , X>0 
1, X < 



(4.15) 



According to (3.9), with I scaling with N such that I := l/N is fixed, we have a = b = N — I, 
and thus cosip = 0, lim^oo cos 7 = 1 — 1, which together imply 



1 1 ,. „, 1 

km Mi = - 

2 2 i— >oo 2 



1 



;i - ly 



(4.16) 



11 



With I identified as S, this is precisely xq(S) in (4.10). We want to use Proposition 
show that with an appropriate choice of scaled variables as implied by (4.10) and (4.1 



4.1 to 



), the 



summations (3.13) defining the correlation kernel K are slowing varying in the summation 
index I and as iV — > oo are given by explicit Riemann integrals. The latter we will recognise 
as specifying the well known Airy process from random matrix theory [151 [11] . 111 addition, for 
purposes of application to the scaled limit of the probability E we need to demonstrate that 
the error term is integrable to the right of the soft edge. 



Proposition 4.2 Let 



x = x (S x ) + XoN- 2 '\ y = xo(S y ) + YaN^ 3 , 
S X = S + stN~ 1/3 , S y = S + trN~ 1/3 , 



with S 7^ 1, xq(S) as specified below (4-10) and a, r such that 



a" 



(i-sy- 



16y/S(2-S) 



T 



(1-S) 2 y/S(2-S) 



2a 



With K as in (3.13), define K by 

K(X,s;Y,t) = Jim aN~ 2 / 3 K(x, NS X ; y, NS y ). 

We then have 

F(X,s 



N->oo 



K(X,s;Y,t) 



i 'v^T / e (i " s)u Ai (u + X) Ai (u + Y)du, s>t 
- F ^ S ) [ e (i_s) "Ai (u + X) Ai (u + Y) du, s<t 

_. F(X, g) AiryP 
' F(Y,t) ^,*,r,tj 



and consequently 



P(fc)({(X,-,Si)}i=i....,fc) := lim (o-N-^p^dx^NS^j^k) 

iv— >-oo 

= det[K AiryP (X J , Sj ; X h s,)]i,J=i,..,A- 



Moreover, 



lim (aN' 

N-+oc 



-2/3 \fc 



dx\ 



in 



dx k P(k)({(xj , NS Xj ) } j=i,...,k) 



Yi 



dXf - / c?A fc p (A . ) ({(A i ,s j )} j= i i ... 5 A ; ). 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



(4.22) 



Proof. Let us introduce l x according to S x = l x /N. Substituting x as given in (4.17) 
into the uniform asymptotic expansion for the Jacobi polynomial (4.15), with n = l x — j, 
a = b = N — l x , gives 



p(N-l x ,N-l x ), x 



jy(N-l x ,N-l x ) 



2x N - l *+ 1 (l 



X 



\N-l x + 



-lM(X*) + 0(N~^)l e ~ XV2 > ^>0 



(4.23) 
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where, in the notation of (4.13), X* is given by 

2 



X* 



Mi _,■ - 1 



■x-3 



(4.24) 



Recalling the definition (3.13) and using the expansion (4.23), K(x,l x ;y,l y ) is given by a sum- 
mation over j of 



jy{N-l x ,N-l x ) 



(-1)1.-1. j ( X (l - X ))N-l,-l ^ 

( y (l _ y ))iV-^+l — j')! A tfN-W-ly) V K l*-j K h-3 a l*-j a h-j 



(4.25) 



x 



( Ai(X*)Ai(F*) + 0(N' 



■2/3\ 



e -x*/2^ x* > 
1, X* < 



i, r* < o 



In fact the summand (|4.25|) is a slowly varying function of j. To see this we first note from 

3 (S x -l) 2 



(4.13) that 



Mi 



lx-j 



— \l S x (2 — S x ) — 



NS X (2-S X 



+ m&x(0{N- 1 ),0{f/N 2 )^ (4.26) 



and v\ x -j = 2oN 2 / 3 + 0(iV x ). This substituted in (4.24) suggests we introduce a continuous 
summation label w by j = wN 1 / 3 to obtain the large N form 



x* = -w { s 1)2 = + x + cKiv- 1 / 3 ) 

2ay/S(2 - S) 



(4.27) 



Furthermore, in terms of w use of Stirling's formula shows that the first line of (4.25) (the 
factors before the Airy functions) reads 



F(X, s) ( 8N^ 3 a 



F(Y, t) VC5 1 - I) 2 .' (XP 
with the error term uniform in w and where 



(s — t)rw 
5(2 - S) 



(1 + 0(A^ 1/3 )) 



(4.28) 



F(x, s ) = (-iy*[x(i-x)r-^N^ {S (2 - g))^ 2/3 ^ e 4gy^ /3 + s3T3(2 : 2S+52) 3^( 2 ^) 2 

(4.29) 

Thus, to leading order the summand depends only on the continuous summation label w, and 
thus can be converted to a Reimann integral to give 



aN- 2 / 3 K{x,NS x ;y,NS y ) 



F(X, s) ( 8a 2 



exp 



(s — t)rw 
S(2 - S) 



x Ai —w- 



(5-1) 



2a^/S(2-S) 
+ 0(N-^ 3 )0(e^ x+Y ^ 2 ). 



+ X Ai -in 



(5-1) 



2^5(2 - S) 



A change of variables 



u 



- w (s-iy 



2a^/S(2-S) 
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+ Yj dw 
(4.30) 

(4.31) 



turns this into 



lim aN- 2/3 K(x,NS x ;y,NS y ) 

2V— >oo 



F(X,s) 
F(Y,t) 



exp 



(t — s)t2<ju 



^s(2-s)(i-sy 

x Ai (u + X) Ai (u + Y) du 



(4.32) 



and the choice of r as in (4.18) completes the proof in the case s >t. Furthermore, essentially 



the same working holds true for the remaining case s < t. 



With the form of K now established, the formula (4.21) for the correlation function follows 



by noting that upon substitution in (3.12) the factor F(Xj,Sj)/F(Xi,si) cancels from the 
determinant. Finally, the validity of (4.22) follows from (3.12) combined with the error term 
in ((4301). □ 



We can make use of Proposition 4.2 to compute the scaled limit of the probability E . Thus 
we know from [18], [3] that the convergence of the integrals (4.22) implies that the scaled limit 



can be applied term-by-term in (3.15). Since the limiting correlations are determinantal, the 



resulting sum can again be summed by appealing to Fredholm integral operator theory (see 
e.g. §9.1]). 



Corollary 4.3 Let Vj and Vj be related as implied by the first line in [J^.ll), and Sj and Sj be 
related as implied by the second line. We have 

lim E ({NS f , (v u l)} i= i,...,„) = det (l - if AiryP ({ Sl , {V h oo)} i=1 ,..., n )) (4.33) 
where K AiryF ({si, (Vi, oo)} i=1) . „ )fl ) is the n x n matrix Fredholm integal operator with kernel 



ir Ai ^ p (X,F;{ S ,;(^,oo)} l=1 ,... in ) = Xxeiv^K^ (X, s f , Y, s 3 ) X y 



i,j=l,...,n 



(4.34) 



The expression (4.33) is precisely that for the cumulative distribution function of the scaled 



largest eigenvalue in the Dyson Brownian motion model of complex Hermitian matrices [16J 
(see Ch. 11] for an account of the model). We remark that the statistical system defined 
by the scaling of the eigenvalues in the Dyson Brownian motion model of complex Hermitian 
matrices about the largest eigenvalue is referred to as the Airy process [TB] . And as noted in the 
Introduction, the Airy process also occurs in random matrix theory as the correlation functions 
for the GUE minor process at the soft edge [TO] . 



5 Universality 

We know turn our attention to the case of a general absolutely continuous distribution T> for the 
height variables, characterized by the corresponding PDF h(x). We have already seen that the 



change of variables (2.5) maps the joint PDF for configurations of interlaced heights sampled 
from T> to the joint PDF for configurations of interlaced heights sampled from the particular 
case T> = E/[o,i] - - the uniform distribution on [0,1]. And in the previous section we have 
identified as the Airy process a particular scaling limit of the joint distribution of the heights 
along y — 1 in the case T> = i7[o,i]- Our aim in this subsection is to show that the change of 
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variables (2.5) implies that there is a choice of scaled variables which also gives the Airy process 
as the scaling limit of the joint distribution of the heights along y — 1. 

Let N be large, and scale the vertical lines in the corresponding bead process so they 
are labelled by a continuous parameter S as specified at the beginning of §4.1. Again in the 

It follows 



bead process picture, (4.1) gives the support of the density in the case V = U\ 



[o,i]- 



immediately from the change of variables (2.5) that the along line S the interval of support for 



a general absolutely continuous distribution with PDF h is specified by the equations 

1 1 rcs(h) i i rd s (h) 

- - -y/S(2-S) = J h(t) dt, - + -y/S(2-S) = J h(t) dt. 



(5.1) 



We note that ds(h) must be finite for S ^ 1. Recalling (4.10) in the case V = £/[o,i] for a scaling 
of positions in the bead model/ heights in the directed solid-on-solid model about the upper 
edge of support, we introduce the renormalized parameter a such that 



d s (l) + XjaN- 2 / 3 = [ 

Jo 



ds^+Xjd-N- 2 / 3 



h(t) dt. 



(5.2) 



For h(t) continuous we ca expand the RHS to leading order N 2//3 and so deduce that to leading 
order in N, 

a = ah(d s (h)). (5.3) 

Hence we conclude that for h(t) continuous and non-zero at t — ds(h), the scaling of positions 
about the upper edge of the support 



x 



dsM + Xj/WsM)^ 3 ) 



(5.4) 



is consistent with the change of variables (2.5) required for the equality of joint probabilities 



(2.6). Thus with the change of variables (5.4) the scaling limit of the joint distribution of 
heights in the general case must be the same as in the particular case V = t/[o,i], so giving our 
sought universality result. 



Corollary 5.1 Let Vj and Vj be related as for Xj and Xj in (5.4)- Let Sj and Sj be related 
as given in the second line of (4.11), where again it is required S ^ 1. For the directed solid- 



on-solid model with heights sampled from a general absolutely continuous distribution T> with 
corresponding PDF h(x), h(x) itself continuous and assumed non-zero at the boundary reference 
point x = ds(h), the limit formula (4-33) remains valid. 



The requirement that S ^ 1 is crucial (S is the reference line in the bead process picture 
about which, on separations of order iV 2 / 3 , the correlations between heights along y = 1 are 
being measured). To see this, let us consider the fluctuation of the height xn,n in the top right 
corner of our directed solid-on-solid model. In the bead process picture, corresponds to the 
particle with the largest coordinate on the line 5 = 1. Now according to the definition of the 
model, xn,n is the largest of the heights. Hence, for a:jv,jv to be less than X we must have that 
all the iV(iV + l)/2 height variables must be less that X, and so 



Pr(x N<N < X) 



x 



h(t) dt 



N(N+l)/2 



(5.5) 
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(here the interlacing requirement plays no role). 
In the case T> = LVii we see from this that 



lim Pr(x NN < 1-X/N 2 

N^oo 



-X/2 



(5.6) 



On the other hand, for T> the exponential distribution it follows from (5.5) that 



lim Vi(x N N < X + 2 log N) = e 

N^oo ' 



(5.7) 



Thus for the fluctuations of the maximum height, both the scale required to get a well defined 
limiting distribution, and the limiting distribution itself, are dependent on T>, in contrast to 



the situation exhibited in Corollary 5.1 
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